d-electron induced icosahedral growth in strontium clusters 
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Strontium clusters of size 2 to 20 atoms have been studied within the local density approximation 
using the ab-initio molecular dynamics method with a plane wave basis in which the s and p 
electron wavefunctions are well represented. The resulting low lying isomers are optimized using a 
full-potential linear mufhn tin orbital method which treats the d-electrons also on an equal footing. 
We find that the d electrons contribute significantly to the binding energy as the cluster size grows, 
expedite the onset of metallicity and induce icosahedral growth. These results are in complete 
agreement with the experimental finding of rare-gas like magic clusters of group IIA elements and 
suggest that the observation of icosahedral growth in some transition metal clusters may have a 
similar origin. 



PACS numbers: 36.40. Cq; 36.40.-c; 71.15.Pd 

An interesting question which has intrigued researchers 
for a long time is how atoms behave collectively as one 
goes from molecular state to the bulk. Small clusters of 
a variety of materials exhibit quite different properties as 
compared to the bulk and are, therefore, attracting much 
attention. It is indeed a challenge of cluster science to 
understand the evolution of the atomic structure and the 
growth of the physical and chemical properties of matter, 
from embrionic state to the bulk. Clusters of group IIA 
elements exhibit a non-metal - metal transition from a 
weakly bonded dimer to a strongly bonded bulk as the 
size grows. Theoretical studies on Be |Q] as well as Mg 
^ clusters showed their magic numbers to be different 
from those of rare gases. These could be qualitatively 
understood in terms of the filling of the electronic shells 
in a spherical jellium model, where one assumes that the 
valence electrons of the atoms move in a jellium of posi- 
tively charged ions [||. Magic clusters of rare gases with 
13, 19, 23, 26, 29, 32, 34, ... atoms, however, correspond 
to icosahedral growth which has been understood to 
arise due to a spherically symmetric pairwise interaction 
and corresponds to a close packing of hard spheres. But 
as one goes further down the group IIA, one encounters 
elements Ca, Sr and Ba where the unoccupied c?-states 
could play an important role in the growth of the atomic 
and electronic properties of their clusters. 

Currently there is little understanding of the growth 
and other properties of metal clusters involving d elec- 
trons, in spite of their importance in several technologi- 
cal applications. Experiments on clusters of Ni ^ have 
provided indirect evidence of icosahedral growth akin to 
rare gas clusters, even though bonding in the latter case 
is very different. Also the abundance spectra of Eu 
Yb [0, Ba ||] and Sr ^ clusters show rare gas type 
magic behavior. However, the origin of this seemingly 
icosahedral growth is not understood. All of these el- 
ements have a ns'^ closed shell electronic configuration, 
and a low-lying empty valence (n — l)d shell. Therefore, 
their clusters can be expected to exhibit similar behav- 
ior. More interestingly, clusters of Eu have been found 



to show different abundances depending upon the nu- 
cleation conditions j^]. Low temperature studies on Eu 
clusters show rare gas type magic clusters whereas exper- 
iments with higher temperature of the carrier gas led to 
abundances expected from metal clusters with a free elec- 
tron like behavior Also, clusters of coinage metals are 
found to exhibit features which appear to be similar to 
those obtained for alkali metal clusters, but deviations 
have also been noted which are believed to suggest an 
icosahedral growth and could be due to the presence of 
d electrons 0,|l|- 

In this letter, we show that the d-electrons induce 
icosahedral growth in metal clusters and explain in par- 
ticular the magic behavior observed in the abundance 
spectra of Sr and Ba clusters. The abundance spectrum 
of small Sr clusters [|| is similar to the one of rare gases 
but it shows Sr II to be magic which can not be under- 
stood from either the jellium model or the rare gas like 
picture. Further, Sris is surprisingly not magic. The 
abundance spectrum of Ba clusters is similar but Bai^ 
is also magic. To understand these, we have carried out 
detailed investigation on Sr clusters. 

We have used the Car-Parrinello (CP) jlj] method and 
simulated annealing technique to obtain the low lying iso- 
mers within the local density approximation. As a first 
approximation to the electronic structure of Sr clusters, 
we ignored the d states in order to save the computer 
time because we expected their occupancy to be small 
in clusters. The ionic potential was represented by the 
norm-conserving pseudopotential of Bachelet et al IT^ ] 
in separable form and the wave functions were expanded 
in a plane wave basis with a cut-off of 8 Ry. The sim- 
ulations were carried out using the s-nonlocality, which 
treats only the s- and p- ionic pseudopotentials accu- 
rately. The cluster was placed in an fee cell of 64 a.u. 
with periodic boundary conditions. This size of the cell 
was large enough so that interaction between the clus- 
ter and its periodic images was negligible. The Brillouin 
zone was sampled by the T point and all calculations 
were performed for the singlet state. For clusters with N 
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< 7, we, however, relaxed a few selected geometries with 
steepest descent method. Structures of the low lying iso- 
mers were subsequently optimized using the full potential 
linear muffin-tin orbital (FPLMTO) method Q which 
treats the s-, p- and d-states on an equal footing. This 
procedure thus allowed us to study the role of d-orbitals 
in the bonding of Sr clusters. 




Figure 1. Minimum energy structures of Sr clusters cal- 
culated using the CP method. 

The lowest energy structures obtained from the CP 
method are shown in Fig. 1. For N < 11, these are sim- 
ilar to the results obtained earlier for magnesium clus- 
ters 1^ except for Srg, which is an octahedron. S'rg is 
a capped pentagonal bipyramid and it is lower in energy 
than a bicapped prism as well as a hexagonal bipyra- 
mid. For Srg, a tricapped prism is lower in energy than 



a double-capped pentagonal bipyramid. Srio is a tetra- 
capped prism whereas in the case of Srn all the five faces 
of a prism get capped and this is lower in energy than an 
icosahedron with two atoms missing. 

For Sri2 simulated annealing yields a tetrahedrally 
packed structure to be of lower energy than an icosa- 
hedron with one vertex atom missing. However, for Sri^ 
a regular icosahedron is of lowest energy, though neither 
Bei3, nor Mgi^ has an icosahedral structure. This is 
due to the fact that the 1/ state (of the jellium model) 
splits into a 3-fold and a 4-fold state. The 3-fold state is 
completely occupied and it is separated from the unoc- 
cupied 4-fold state by a gap of about 0.5 eV. The lowest 
energy structure for Sri4 is a capped icosahedron. Sri^ 
is an incomplete icosahedron with 3 caps while S'rig is 
a distorted icosahedron with 3 capping atoms. Srn is 
a highly symmetric Frank-Kasper polyhedron with 16 
atoms forming a cage and one atom at the center. For 
Sris, the structure has similarity with Srii and for Sr ig 
an ftcp-like structure (see Fig. 1) has lower energy than a 
double-icosahedron. Sr2Q is a double icosahedron with a 
cap in the middle. These results clearly show that there 
is no systematic icosahedral growth if the system is s — p 
bonded. 
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Figure 2. (a ) Binding energy per atom and (b ) second 
derivative of energy for Sr clusters as a function of clus- 
ter size. Solid (broken) lines denote the results from 
FPLMTO (CP) calculation. 



2 



The binding energy per atom is shown in Fig. 2. It 
is noted that the approach to the bulk value (« 2.0 
eV/atom is fairly weak. The second derivative of 
energy, AEn = E^+i + i^Ar-i — 2i?jv, Em being the en- 
ergy of a TV atom cluster, is significantly positive for N 
— 4, 7, 10, 15 and 18 favoring a magic behavior for these 
clusters (see Fig. 2). An interesting result to note here 
is that although the lowest energy structure for Sri^ is 
an icosahedron, it is not magic. On the other hand, clus- 
ters with 10 and 18 atoms are not magic experimentally. 
Thus these results do not agree with experiments on Sr 
clusters. 

In order to investigate the role of c?-electrons, the low- 
lying isomers of all the clusters were optimized using the 
FPLMTO method, with s, p and d basis orbitals. The 
resulting lowest energy structures are shown in Fig. 3 
for cases in which the structures are different from those 
obtained from the CP method. It is found that the d elec- 
trons lead to a reduction in the bond lengths in all the 
structures investigated. Most of the larger clusters which 
did not have an icosahedral structure to be of lowest en- 
ergy in the CP method, now turned out to have very 
well-defined icosahedral geometries. Other structures re- 
tained the geometry but there is relaxation in the bond 
lengths which have a spread in their values. S'rg now 
turned out to be a pentagonal bipyramid with two caps. 
Similarly, Stiq, Sru and Sr 12 now yielded pentagonal 
bipyramid with three, four and five caps respectively. 
These can also be visualized as icosahedra with three, 
two and one atom missing. From Sri^ onwards all clus- 
ters are capped icosahedra until one obtains a double- 
icosahedron for Srig with the exception of Sri<^ which 
is a bicapped hexagonal antiprism (see Fig. 3). How- 
ever, the bicapped icosahedron, at the binding energy of 
1.591 eV/atom, is nearly degenerate with the minimum 
energy structure at 1.595 eV/atom. It is clear that now 
the icosahedral growth becomes more favorable. Similar 
structures have been observed in the experimental stud- 
ies on Ni clusters , and also in theoretical studies using 
empirical many-body potential for Ni atoms [^|. 

The binding energy is substantially larger when the 
d electrons are also included (Fig. 2). For Sr2 it is 
0.29 eV/atom in fair agreement with the value of 0.27 
eV/atom obtained by Ortiz and Ballone jl^. However, 
for SriQ the binding energy is 1.75 eV/atom. Therefore, 
we find a transition from a weak bonding to a strong 
metallic bonding as a function of the cluster size. It 
is also interesting to note that the contribution due to 
the inclusion of the d-states to the binding energy per 
atom is very small (about 0.18 eV) for a dimer, but it 
grows to approximately 1.1 eV for the 19-atom cluster 
I p^ . Thus, the c?-states become increasingly important 
as the sp — d hybridization increases with size. The pres- 
ence of c?-states expedites the approach to the bulk prop- 
erties as also evidenced from Fig. 4 by a small highest 
occupied (HOMO) and the lowest unoccupied molecular 



orbital (LUMO) gap (« 0.1 eV for N > 10). 




Figure 3. Lowest energy structures of Sr clusters ob- 
tained from FPLMTO method, which differ from those 
obtained by the CP method. The d-states favor an icosa- 
hedral growth governed by capping the Sr^ pentagonal 
bipyramid. 

We consider our results as compelling evidence that the 
presence of unoccupied c?-states make icosahedral pack- 
ing more favorable in Sr clusters. This is the closest 
packing of atoms in clusters. Thus the reason for the oc- 
currence of rare-gas like magic numbers in the abundance 
spectrum of Sr clusters is the presence of unoccupied d 
states close to the s level which participate very signifi- 
cantly in the bonding. 

The second derivative of energy now shows that Stiq 
is no longer a magic cluster. On the other hand, Sru 
as well as Srig now become magic in agreement with the 
experiment. This substantiates the generally observed 
feature that the s-p bonded metals fall into the jellium 
picture. It shows in addition that the presence of c?-states 
renders the jellium picture inappropriate. Our results 
with c?-electrons also show Sri^ not to be magic which is 
also in agreement with the experimental data. We find 
that for Sr II there is a significant contraction in the bond 
length which leads to a stronger bonding contribution 
from d states. Whereas, for Sri^, the bond lengths are 
slightly longer as compared to the one in Sru. 

Thus our results agree with the experimental finding 
for Sr and suggest that the observation of rare gas like 
magic behavior for Ba, Yb and Eu are also due to the sig- 
nificant contribution of d-states to the bonding in these 
clusters. Since 5d states get more populated in Ba, Yb 
and Eu, we would expect these to make a larger con- 
tribution which may give rise to a magic behavior for 
13-atom cluster also as is found in the case for Ni clus- 
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ters. Our results also suggest that the observation of 
temperature dependent abundance spectra could also be 
due to a change in the sp — d hybridization resulting from 
thermal expansion in clusters. Note that for weak contri- 
bution from d-electrons we find 10- and 18-atom clusters 
to be magic which agrees with the experiments on Eu. 
Indeed, in our CP simulations on Sr clusters we observed 
large increase in the interatomic distance at finite tem- 
peratures which should make such clusters susceptible 
to structural changes. We feel that the icosahedral type 
of growth observed in coinage metal clusters be- 
yond a size of 40 atoms may be due to a growing sp-d 
hybridization. It is also appropriate to recall here the 
structural changes in gold clusters due to electron irra- 
diation ||l9|, which may be a heating effect of radiation. 
It is known that the s — d hybridization plays a crucial 
role in the stability of the fee phase of gold. Therefore, 
there is a critical balance arising from d electrons which 
transform the clusters to a closed packed structure. 
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Figure 4. The HOMO-LUMO gap for Sr clusters. One 
can see that it is not correlated with the peaks of AE^. 

To summarize, we have calculated the structure and 
electronic properties of Sr clusters in the range of 2-20 
atoms, using the CP and the FPLMTO methods. We 
find that the d-electrons contribute substantially to the 
binding energy and induce icosahedral growth in these 
clusters. The contribution increases from a very small 
value for the dimer to a large one for a 20-atom cluster. 
The calculated magic numbers are in agreement with the 
experimental results. These are the first results where 
the effect of d electrons has been studied on the growth 
of clusters and we hope that our studies would stimulate 
further work in the understanding of the properties of 
transition metals clusters. 
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